Supersymmetric quenched complexity in the Sherrington-Kirkpatrick model 
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o : 

^ By using the BRST supersymmetry we compute the quenched complexity of the TAP states in 

the SK model. We prove that the BRST complexity is equal to the Legendre transform of the static 
free energy with respect to the largest replica symmetry breaking point of its overlap matrix. 
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. A key issue in the physics of complex systems is the computation of the entropy of the metastable states, normally 

' called complexity in spin-glasses, and configurational entropy in structural glasses and supercooled liquids. A knowl- 
edge of the complexity is crucial for understanding the dynamics of a system when this is heavily influenced by strong 
metastability effects. Moreover, in some theoretical frameworks, the drop in the number of accessible states leads to 
^ ■ an ergodicity breaking transition. In this context the complexity is essential also from the thermodynamic point of 
D ■ view, as in the Adam-Gibbs theory of the thermodynamic glass transition 0. 

Despite its enormous theoretical relevance, there are few analytic calculations of the complexity in glassy systems, 
and this for a very good reason. In a nutshell, to find the complexity we have to compute the number of local 
minima (metastable states) of some state function, which is typically highly nontrivial. Just to fix ideas, we may 
think that this function is the Hamiltonian H. To compute the complexity, we must impose that the gradient of H 
(the force) vanishes in the local minima, and we have to include as a normalization factor the second derivative of 
H (the Hessian). Moreover, we may want to classify the metastable states according to the value of H, that is to 
their height in the landscape. Therefore, beside the force and the Hessian, we must include the state function itself 
in the calculation. Computing the complexity is thus a formidable technical task, since we have to deal with three 
fH , very complicated functions: H , dH, and d^H. In comparison, the calculation of the partition function, which just 
Q involves H, is an easy business. 

O , This apparent difficulty in the calculation of the complexity stems from the fact that most methods treat H, dH 
and d^H as three independent objects, when of course they are not. Every calculation which fails to capture the fact 
—i ,, that it is essentially just one function, iJ, that we are dealing with, effectively wastes a crucial information. It would 
^ ■ be therefore important to find a tool which exploits this information to simplify the calculation of the complexity. 
The Becchi-Rouet-Stora-Tyutin (BRST) supersymmetry 0,01 seems to be such a tool. As it was first noted in for 
a particular model, a BRST calculation of the complexity is in fact equivalent to the one of the partition function. 
, This is indeed what we expect from a method which does not treat H, dH and d^H as independent functions. 
' The formal equivalence between complexity and standard thermodynamics found in 0| by means of the BRST 
. supersymmetry is a very important theoretical issue. In the context of spin-glasses the existence of such a connection 
■ has been much investigated in the past 00000 0, O, 0, ^3 ■ In a classic paper , Bray and Moore first 
calculated the complexity of the Sherrington-Kirk patr ick (SK) model by counting the number of local minima of 
the Thouless- Anderson-Palmer (TAP) free energy which in mean- field spin-glasses is the state function discussed 
above. The same authors also noted in some deep formal connections between TAP complexity and standard 
'"^ ■ thermodynamics, while De Dominicis and Young showed in that TAP and static approaches were in fact equivalent, 
, once some key hypothesis were made. These studies culminated in a remarkable work where Bray, Moore and 
O ■ Young uncovered a sort of Legendre transform relationship between TAP complexity and static free energy. 

A method to compute the complexity which does not rely on the existence of a TAP free energy, was introduced 
by Monasson and by Franz and Parisi The basic idea is to introduce a coupling between different systems 
forcing them to live in the same metastable state. The free energy cost of such a constrained super-system is equal 
^ • to the entropic contribution of the metastable states, that is the complexity. Within this approach, close connections 
9^ between complexity and thermodynamics, similar to those found in the TAP context in 0, were found. In particular, 
in spin-glass models with one step of replica symmetry breaking [l7l | , the formulation of Monasson shows that the 
complexity is equal to the Legendre transform of the static free energy with respect to the breaking point x of the 
overlap matrix ID. l.'i . 

Despite all these investigations, it is fair to say that a general formal connection between complexity of the 
metastable states and static free energy, has not been proved yet. In particular, it is unclear how the Legendre 
transform method of should be used in systems with more than one step of replica symmetry breaking, as the SK 
model. In fact, none of the previous SK investigations 00 succeeded in proving the existence of a sharp Legendre 
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transform relationship as in one-RSB systems. 

In this Letter we find for the first time an exact connection between complexity of the metastable states and 
static free energy in the SK model: we prove that the quenched TAP complexity obtained by means of the BRST 
supersymmetry is the Legendre transform of the static free energy with respect to the largest breaking point of its 
overlap matrix. Our result confirms the validity of the Legendre transform method of |lOl IllL Il2| , and its consistency 
with the investigations of [ESS Hi- Moreover, our findings strongly suggest that the BRST supersymmetry should 
be considered as an essential tool for the computation of the complexity in more general glassy systems. 

The complexity of the TAP states with free energy density /, at inverse temperature /3, is defined as 0, 

1 ^ ^ r ^ 

/) = ^ log 5^ <5 [NPf - l3FTAp{m^)] = JJ^og dr e^''^^ ^ g-r/JF^^H™") , (i) 

where m" = {mf}, are the local magnetizations at site i = 1 . . . A^, in state a — 1 . . .N. A state m" is defined as a 
local minimum of the TAP free energy Ftap{^) 15|. If we define the thermodynamic potential r), 

M 

exp {-l3Nr *) = ^ g-r/3F^^^(m°) ^ ^2) 

a=l 

we can use the steepest descent method in (^^l, and obtain the complexity as the Legendre transform of 5'(/3,r), 

S(/3,/)=/3r/-/3rvI/(/3,r) , (3) 
where the parameter r = r(/3, /) is fixed by the equation, 

^^P^r)+r^^^^=f . (4) 

From (jSJ we see that for r = 1 the potential 5" must be equal to the standard static free energy of the system F{P), 
calculated in the TAP context. This calculation was first performed in where it was shown that the relation 
^'(/3,r = 1) = F{l3) only held if some suitable assumptions were made. In it was proved that the assumptions 
used in Q were in fact a general consequence of the BRST supersymmetry. However, the situation was less clear for 
generic values of r, since the calculations of 0, B] for r 7^ 1 seemed to explicitly break the BRST invariance 
what follows we perform a supersymmetric quenched calculation of r), and prove that this potential is intimately 
related to the static free energy F{(3) even for r 1. The TAP free energy for the SK model is given by [l^ . 

f3 1 

f3FTAp{m) = - 2 X] J^j^i^o + ^^Min-i) , (5) 



with, 



1 _ 
(j}o{m) ^ - log(l - m^) + m tanh"^(TO) - log2 - -^(1 - . (6) 

The variable q is the self-overlap of the TAP states, q = X^i The quenched couplings J are random variables 
with Gaussian distribution and variance N. From |2Jl we have that the quenched potential \E'(/3, r) is, 



-/?r«'(/3,r) = — logp(/3,r|J)" , (7) 
nN 



with, 

AT 



p(/3,r| J) = e-rPF^AAm-) ^ fY[dm, 5{d,FTAp{ni)) \ Aei{d,d,FTAp{ni))\ ^-PrF^.A^) (g) 

a=l '' i 

In {Tj) we have N ^ 00 and n ^ 0, and the over-bar indicates an average over the disorder. As usual, the modulus 
of the determinant will be dropped. This amounts to assume that at sufficiently low temperatures the largest part of 
TAP solutions are minima 18]. Of course, any method which drops the modulus is doomed to fail if stable minima 
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are subdominant with respect to unstable saddles. After introducing the commuting fields Xi to implement the delta 
functions, and the anti-commuting (Grassmann) fields V'iiV'j for the determinant, we find. 



where the action S is given by, 

S{m,x,ip,'ip) = ^XidiFTAp{m) + ^p^ipjdidjFTApim) - rFTAp{m) 

i ij 

By averaging p{P, rj J)" over the disorder we obtain the following effective action, 
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where (pi — dm4>o and (j)2 — d^^ipQ. The scalar overlap q will now be generalized by introducing the overlap matrix 
Qab = m'^ • TO^- This form of the action is different from the one in where the delta function enforcing the TAP 

equations was used to eliminate the factor J^m^m" in pU|l . This is a crucial point: this method explicitly breaks the 
BRST invariance of the action, which is the crucial tool to establish the exact connection with the static free energy. 
We therefore do not use this method, and keep the whole BRST invariant action p^ . 



In order to linearize the quadratic terms we introduce the usual Lagrange multipliers, m^m* — > Xab, m'^x'' 



Wab, 



^a^b After this is done, the integrals in x and ip,ip become Gaussian and can be performed explicitly. 

Moreover, the action factorizes and for ^ cxo we can use the steepest descent method, to get. 



-Br "if (13, r) = lim - 

n-»0 n 



Eo + log y" n 



Following 0, 0, we define, 



Bab = f3Hl-qaa)Sab + t''' 
Aab = -/?'(l-gaa)<5ab-U;°' 



(13) 
(14) 



and therefore obtain (details will be given elsewhere), 
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qi, + X^'qab ~ rA^'q^b 



£{m-) = -r V MQaa, m'^) + V A'^^m'^m" + log det ( + B, 



■log[(2^/32)"detg,b 



(15) 
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(16) 



The parameters Aab, Bab, ^ab and qab must be fixed by the saddle point equations, and it is easy to show that Bab — 
is solution. It is important at this point to consider the role of the supersymmetry. In 4] it was noted that action 
()10|l is invariant under a generalization of the BRST supersymmetry j^, |3 , 
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where e is an infinitesimal Grassmann parameter. If we calculate the variation of niiipi and Xiipi [3. Il8|. we obtain 
the two BRST equations, 



After some algebra, these equations become, 



(V'iV'i) + {rriiXi) = 
r {"ipiipi) + {xiXi) = 



1 



(17) 
(18) 



(19) 
(20) 



and it is possible to show that the remaining saddle point equations are indeed satisfied by H19|) and (|20|) . The only 
saddle point equation we are left with is obtained by doing the variations of ifTs)) and l(TB|) with respect to Xab- This 
gives. 



(21) 



where the average is performed with the distribution exp(£(TO")). If we use equations (|19|1 . H20|l and H21|l into H15|l 
and Hl()ll . and make the change of variable to" h°- = tanh^^(TO''), we finally obtain. 
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Expression (|22|l is different from the one of 0, Isl : the BRST supersymmetry automatically selects one saddle point 
in the space of parameters, and in so doing it drastically reduces the number of parameters, compared to 0,13. The 
computation of '^{j3,r) has at this point the same degree of difficulty as the one of the standard free energy F{f]), 
with just one overlap matrix qab to be fixed variationally. We shall now show that the connections between Vl/(/3, r) 
and F((3) are in fact much deeper than that. 

The general form of the quenched free energy in the SK model is [l^ . 
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where Qaf3 is the ris x ns overlap matrix, with Ug 0. If in (|22|l we use the relation, 
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we obtain, 
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The trace terms in equations (|23|) and (|25|l suggest the relation ns = r ■ n. Once this identification is done, we can 
connect the u" spin variables {a = 1, . . . , n^), to the spin variables (a = 1, . . . , n; = 1, . . . , r) in the following 
way, 



(cTi, . . . , cr„ 
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Let us now assume that the potential 'I'(/3,r) (and thus the TAP complexity) is calculated at k levels of replica 
symmetry breaking (RSB) 17] . The TAP overlap matrix qab is then given by, 



fc+i 



(26) 
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with Uk+i = 1 and £^b'^^ = ^ab- The matrices e("'i'») are nxn ultrametric block matrices, equal to one on the diagonal 
blocks of size j/i, and zero elsewhere. The variables yi are thus the replica symmetry breaking points. In the TAP 
approach the diagonal of the overlap matrix, qaa = 9fc+i, contains the self-overlap of the states, and for this reason 
Uk+i = 1. There are k + 1 values of the overlap, but only k nontrivial breaking points, and thus qab is a A:-RSB matrix. 
Given this form of qab, it is possible to prove (details will be given elsewhere) that, 

n r n r rn 

E E -EE - E qT -"-/^ ' (27) 

ab fiv a fj. 

where Q)^q is a standard rn x rn RSB matrix, with k + 1 levels of replica symmetry breaking. More precisely. 



fc+i 

QiT =10 + E(«^ - 9-1) 4"""^ ■ (28) 
1=1 

(k) 

From this formula we see that the entries of are the same as g^^^ , whereas the fc + 1 replica symmetry breaking 

points Xi of Qg^f^ are rescaled by a factor r, that is Xi = r yi. In particular, the largest breaking point x = Xk+i of 
the static matrix Q^'^^^^ is given by, 

x = r . (29) 

By inserting relation H27(l into H25(l . we finally obtain, 

*(/3,H9i'^)=n/3|Qir^) , (30) 

We have thus proved that the thermodynamic potential ^'(/3,r) calculated at the k RSB level is equal to static free 
energy F{f3) calculated at the k + 1 RSB level. The replica symmetry breaking points of the static matrix Q^fs are 
simply the ones of the TAP matrix qab rescaled by the parameter r, and the extra k + 1-th breaking point of Qap 
is equal to r. This rescaling was first noted in 0, and later in although the lack of BRST symmetry of those 
calculations prevented to prove equation (EUJ). 

From equation Q, and given the relation between ^'(/?, r) and F{(3), we finally have the general Legendre equation 
connecting the quenched complexity of the TAP states to the standard static free energy in the SK model, 

S(/3,/)=/3x/-/3xF(/3;x) , (31) 

with the largest breaking point x fixed by the equation, 

/ = mx)+.^^. (32) 

This result can be summarized as follows: the supersymmetric quenched complexity of the TAP states is the Legendre 
transform of the static free energy with respect to the largest breaking point x of its overlap matrix. This result allows 
to compute the exact quenched complexity in the SK model, that is the complexity at the full- RSB level, as the 
Legendre transform of the full- RSB static free energy 20]. 

It has been conjectured in .21] that in systems with more than one step of replica symmetry breaking the complexity 
of clusters at level i is given by the Legendre transform of the free energy with respect to the breaking point Xi. For 
Xi = Xmax clusters are just states, and our result is recovered. It would be interesting to study whether the conjecture 
of |2ll | can be exactly proved within the supersymmetric formalism used here. 

We thank A. Crisanti, L. Lcuzzi, R. Monasson, A. Montanari, F. Ricci-Tersenghi, and T. Rizzo for some interesting 
discussions. 
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